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Abstract 

In this paper, we study the following problem of reconstructing a simple polygon: Given a 
cyclically ordered vertex sequence of an unknown simple polygon P of n vertices and, for each 
vertex v of P, the sequence of angles defined by all the visible vertices of v in P, reconstruct the 
polygon P (up to similarity) . An 0(n 3 log n) time algorithm has been proposed for this problem. 
We present an improved algorithm with running time 0(n 2 ), based on new observations on the 
Oh- geometric structures of the problem. Since the input size (i.e., the total number of input 

visibility angles) is Q(n 2 ) in the worst case, our algorithm is worst-case optimal. 
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1 Introduction 



In this paper, we study the problem of reconstructing a simple polygon P from the visibility angles 
measured at the vertices of P and from the cyclically ordered vertices of P along its boundary. 
q \ Precisely, for an unknown simple polygon P of n vertices, suppose we are given (1) the vertices 

ordered counterclockwise (CCW) along the boundary of P, and (2) for each vertex v of P, the 
angles between any two adjacent rays emanating from v to the vertices of P that are visible to v 
such that these angles are in the CCW order as seen around v, beginning at the CCW neighboring 
| vertex of v on the boundary of P (e.g., see Fig. [T|). A vertex u of P is visible to a vertex v of P if the 

line segment connecting u and v lies entirely in P. The objective of the problem is to reconstruct 
the simple polygon P (up to similarity) that fits all the given angles. We call this problem the 
polygon reconstruction problem from angles, or PRA for short. Figure [2] gives an example. 

The PRA problem has been studied by Disser, Mihalak, and Widmayer [6j, who showed that the 
solution polygon for the input is unique (up to similarity) and gave an 0(n 3 logn) time algorithm 
for reconstructing such a polygon. Using the input, their algorithm first constructs the visibility 
graph G of P and subsequently reconstructs the polygon P. As shown in [6J, once G is known, 
the polygon P can be obtained efficiently (e.g., in 0(n 2 ) time) with the help of the angle data and 
the CCW ordered vertex sequence of P. 

Given a visibility graph G, the problem of determining whether there is a polygon P that 
has G as its visibility graph is called the visibility graph recognition problem, and the problem 
of actually constructing such a polygon P is called the visibility graph reconstruction problem. 
Note that the general visibility graph recognition and reconstruction problems are long-standing 
open problems with only partial results known (e.g., see [1] for a short survey). Everett [8] showed 
that the visibility graph reconstruction problem is in PSPACE, but no better upper bound on the 
complexity of either problem is known. In our problem setting, we have the angle data information 
and the ordered vertex list of P; thus P can be constructed efficiently after knowing G. 
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Figure 1: Illustrating the angle measurement at a vertex v: The angle sequence (01,02,03,04) is given. 



Fi gure 2: The given measured angles are ordered CCW along the polygon boundary (left); reconstructing a simple 
polygon that fits these angles (right). 

Hence, the major part of the algorithm in [6] is dedicated to constructing the visibility graph 
G of P. As indicated in [6], the key difficulty is that the vertices in this problem setting have no 
recognizable labels, e.g., the angle measurement at a vertex v gives angles between visible vertices 
to v but does not identify these visible vertices globally. The authors in [B] also showed that some 
natural greedy approaches do not seem to work. An 0(n 3 log n) time algorithm for constructing 
G is given in [6]. The algorithm, called the triangle witness algorithm, is based on the following 
observation: Suppose we wish to determine whether a vertex Vi is visible to another vertex Vj] 
then Vi is visible to Vj if and only if there is a vertex v\ on the portion of the boundary of P from 
Vi + \ to Vj-i in the CCW order such that v\ is visible to both and Vj and the triangle formed by 
the three vertices V{,Vj, and V\ does not intersect the boundary of P except at these three vertices 
(such a vertex v\ is called a triangle witness vertex). 

In this paper, based on the triangle witness algorithm [6], by exploiting some new geometric 
properties, we give an improved algorithm with a running time of 0(n 2 ). The improvement is 
due to two key observations. First, in the triangle witness algorithm [6j, to determine whether a 
vertex vi is visible to another vertex Vj, the algorithm needs to determine whether there exists a 
triangle witness vertex along the boundary of P from Vi+i to Vj—\ in the CCW order; to this end, 
the algorithm checks every vertex in that boundary portion of P. We observe that it suffices to 
check only one particular vertex in that boundary portion. This removes an 0{n) factor from the 
running time of the triangle witness algorithm [6|. Second, some basic operations in the triangle 
witness algorithm [6] take O(logn) time each; by utilizing certain different data structures, our 
new algorithm can handle each of those basic operations in constant time. This removes another 
O(logra) factor from the running time. Note that since the input size is 0(n 2 ) in the worst case 
(e.g., the total number of all visibility angles), our algorithm is worst-case optimal. 

As shown in [6j, if only the angle measurements are given, i.e., the ordered vertices along the 
boundary of P are unknown, then the information is not sufficient for reconstructing P. In other 
words, it may be possible to compute several simple polygons that are not similar but all fit the 
given measured angles (see [6] for an example). 
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1.1 Related Work 



The problems of reconstructing geometric objects based on measurement data have been studied 
extensively (e.g., [2J [3J [T2], [T31 [15]). As discussed above, the general visibility graph recognition 
and reconstruction problems are in PSPACE [8J and no better complexity upper bound is known 
so far (e.g., see [10]). Yet, some results have been given for certain special polygons. For example, 
Everett and Corneil [9] characterized the visibility graphs of spiral polygons and gave a linear 
time reconstruction algorithm. Choi, Shin, and Chwa [3J, and Colley, Lubiw, and Spinrad [5] 
characterized and recognized the visibility graphs of funnel-shaped polygons. 

By adding extra information, some versions of the problems become more tractable. O'Rourke 
and Streinu [133 considered the vertex- edge visibility graph that includes edge-to-edge visibility 
information. Wismath [T7] introduced the stab graphs which are also an extended visibility struc- 
ture and showed how parallel line segments can be efficiently reconstructed from it. Snoeyink |15] 
proved that a unique simple polygon (up to similarity) can be determined by the interior angles at 
its vertices and the cross-ratios of the diagonals of any given triangulation, where the cross-ratio 
of a diagonal is the product of the ratios of edge lengths for the two adjacent triangles. Jackson 
and Wismath [12] studied the reconstruction of orthogonal polygons from horizontal and verti- 
cal visibility information and gave an 0{n log n) time reconstruction algorithm. Biedl, Durocher, 
and Snoeyink [2] considered the problem of reconstructing the two-dimensional floor plan of a 
polygonal room using different types of scanned data, and proposed several problem models. Si- 
dlesky, Barequet, and Gotsman p3] studied the problem of reconstructing a planar polygon from 
its intersections with a collection of arbitrarily-oriented "cutting" lines. 

Reconstructing a simple polygon from angle data was first considered by Bilo et al. [3], who 
aimed to understand what kinds of sensorial capabilities are sufficient for a robot moving inside an 
unknown polygon to reconstruct the visibility graph of the polygon. It was shown in [3] that if the 
robot is equipped with a compass to measure the angle between any two vertices that are currently 
visible to the robot and also has the ability to know where it came from when moving from vertex 
to vertex, then the visibility graph of the polygon can be uniquely reconstructed. Reconstruction 
and exploration of environments by robots in other problem settings have also been studied (e.g., 
see [2CEI1QI]). 

The rest of this paper is organized as follows. In Section we give the problem definitions in 
detail and introduce some notations and basic observations. To be self-contained, in Section [3j we 
briefly review the triangle witness algorithm given in [6] . We then present our improved algorithm 
in Section 01 

2 Preliminaries 

In this section, we define the PRA problem in detail and introduce the needed notations and 
terminology. For ease of discussion and comparison, some of our notations follow those in [6]. 

Let P be a simple polygon of n vertices Vq, v\, ■ ■ ■ , w n -i i n the CCW order along P's boundary. 
Denote by G = (V, E) the visibility graph of P, where V consists of all vertices of P and for any 
two distinct vertices i?j and Vj, E contains an edge e(i?j,iu) connecting Vi and Vj if and only if 
Vi is visible to Vj inside P. In this paper, the indices of all fj's are taken as congruent modulo 
n, i.e., if n < i + j < 2n — 1, then Vi+j is the same vertex as vi, where I = i + j — n (or 
I = (i + j) mod n); similarly, if —n < i — j < 0, then is the same vertex as vi, where 
I = i — j + n. For each Vj £ V, denote by deg{vi) its degree in the visibility graph G, and denote 
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by vis(vi) = (visi(vi) , vis2(vi) , . . . ,visd eg (vi){ v i)) the sequence of vertices in P visible to V{ from 
v-i + \ to Vi-\ ordered CCW around V{. We refer to vis(vi) as Uj's visibility angle sequence. Note 
that since both Vi-i and Vi + \ are visible to «j, vis\{vi) = Vi + \ and visd eg ( Vi )(vi) = For any 

two vertices Vi,Vj in V, let ch(vi,Vj) denote the sequence (vi,Vi+i, . . . ,Vj) of the vertices ordered 
CCW along the boundary of P from vi to vj. We refer to ch(vi,Vj) as a chain. Let |c/i(t>j,Uj)| 
denote the number of vertices of P in the chain ch(vi,Vj). 

For any two distinct vertices vi,Vj E P, let p(vi,Vj) be the ray emanating from t>j and going 
towards Vj. For any three vertices v,Vi,Vj € P, denote by Z v (vi,Vj) the CCW angle defined by 
rotating p(v,V{) around v to p(v,Vj) {vi or Vj need not be visible to v). Note that the values of 
all angles we use in this paper are in [0, 2n). For any vertex v € P and 1 < i < j < deg(v), let 
Z v (i,j) be Z v (visi(v),visj(v)). 

The PRA problem can then be re-stated as follows: Given a sequence of all vertices vq,vi,..., v n -\ 
of an unknown simple polygon P in the CCW order along P's boundary, and the angles Z v (i, i + 1) 
for each vertex v of P with 1 < i < deg(v), we seek to reconstruct P (up to similarity) to fit all 
the given angles. Without loss of generality, we assume that no three distinct vertices of P are 
collinear. 

It is easy to see that after 0(n 2 ) time preprocessing, for any v G V and any 1 < i < j < deg(v), 
the angle Z v (i,j) can be obtained in constant time. In the following discussion, we assume that 
this preprocessing has already been done. Sometimes we (loosely) say that these angles are given 
as input. 

The algorithm given in [6] does not construct P directly. Instead, the algorithm first computes 
its visibility graph G = (V,E). As mentioned earlier, after knowing G, P can be reconstructed 
efficiently with the help of the angle data and the CCW ordered vertex sequence of P. The 
algorithm for constructing G in [6] is called the triangle witness algorithm, which will be briefly 
reviewed in Section [3J Since V consists of all vertices of P, the problem of constructing G is 
equivalent to constructing its edge set E, i.e., for any two distinct vertices Vi,Vj € P, determine 
whether there is an edge e(vi,Vj) in E connecting and Vj (in other words, determining whether 
Vi is visible to Vj inside P). 

To discuss the involved algorithms, we need one more definition. For any two vertices Vi, Vj € P 
with \ch(vi+i, Vj-i)\ > 1, suppose a vertex vi € ch(vi + \, Vj-i) is visible to both V{ and Vj; then we 
let v'j be the first visible vertex to Vi on the chain ch(vj,Vi) and let v[ be the last visible vertex to 
Vj on the chain ch(vj,Vi) (e.g., see Fig. [3]). Intuitively, imagine that we rotate a ray from p(vi,vi) 
around Vi counterclockwise; then the first vertex on the chain ch(vj,vt) hit by the rotating ray is 
v'j. Similarly, if we rotate a ray from p(vj,vi) around Vj clockwise, then the first vertex on the 
chain ch(vj, Vi) hit by the rotating ray is v[. Note that if Vi is visible to Vj, then v'j is Vj and v[ is 

Vi. We denote by Zl^v^Vj) the angle Z Vi (vi,v'j) and denote by zJ 3 (wj,^) the angle Z Vj (v' i ,vi). It 
should be pointed out that for ease of understanding this paper, the above statement of defining 
is different from that in [S] but they refer to the same angles in the algorithm. The motivation 
for defining will be clear after discussing the following lemma, which has been proved in [B]. 

Lemma 1 ff^ For any two vertices Vi,Vj € P with \ ch{v j+i, Vj-\)\ > 1, Vi is visible to Vj if 
and only if there exists a vertex vi on ch(vi+\,Vj-i) such that v\ is visible to both Vi and Vj and 
Ai(v h Vj) + Aj{vi,vi) + Z Vl (vj,Vi) = IT. 

Since the above lemma is also critical to our improved algorithm in Section 01 we sketch the 
proof of the lemma below. 
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Figure 3: Illustrating the definitions of v[ and v'j, and Zl.(vi,Vj) and ZJ. (vt,vi) which are the angles Z Vi (vi : 
and Z Vj (v[,vi), respectively. 



For any two vertices Vi,Vj G P with \ch(vi+i,Vj-i)\ > 1, if t>j is visible to Vj, then it is not 
difficult to see that there must exist a vertex v\ G ch(vi+i, Vj-x) that is visible to both Vi and Vj. 
Since the three vertices Vi, Vj, and vi are mutually visible to each other, it is clear that the triangle 
formed by these three vertices does not intersect the boundary of P except at the three vertices, 
implying that Zj^Vj, Vj) + Zl^Vi, vi) + Z Vl (vj, vi) = Z Vi (v h Vj) + Z Vj (vi,vi) + Z Vl (vj,vi) = it. Such 
a vertex v\ is called a triangle witness of the edge e(vi,Vj) in E. 

Suppose Vi is not visible to Vj) then it is possible that there does not exist a vertex vi G 
ch{vi + \,Vj-\) which is visible to both Vi and Vj. In fact, if there exists no vertex vi G ch(vi+i,Vj—i) 
that is visible to both vi and Vj, then Vi cannot be visible to Vj. Hence in the following, we assume 
that such a vertex vi exists, i.e., a vertex vi G ch{vi + \, v is visible to both Vi and Vj (yet Vi 
is not visible to Vj). Note that vi is not collinear with Vi and Vj. If Z Vl (vj,vi) > ir (i.e., the 
chain ch(vi,Vj) is part of the boundary of P that blocks the visibility between vi and Vj), then 
the lemma obviously holds. Otherwise, Z v .(vj,Vi) < tt, and in this case, the visibility between Vi 
and Vj is blocked by the chain ch(vj,vi). Since Vi is not visible to Vj, for any choice of such a 
vertex vi, the angle Z Vi (vi,Vj) is not given by vis visibility angle sequence vis(vi). The "closest 
approximation" for Z Vi (vi,Vj) in this case is determined by a vertex v on the chain ch(vj,Vi) such 
that v becomes vj if and only if Vj is visible to v,- L . As in the definition of Z% i (v\ ,Vj), the vertex v'- 
is such a vertex v, i.e., the first visible vertex to Vi on the chain ch(vj,Vi). Similarly, the vertex 
V{ in Z Vj (vi,V[) is "replaced" in the definition of zX j {vi,vi) by v[, i.e., the last visible vertex to vj 
on ch(vj,Vi). Clearly, when Vj is not visible to Vi, it must hold that Z Vi (vi,Vj) < Z Vi (v\,Vj) and 

Z^v'^vi) < Z Vj (vi,vi). Therefore, zt.(vi, Vj) + /l^Vi, v{) + Z Vl (vj, vi) = Z Vi (v h u<) + Z Vj (v' i: vi) + 
Z Vl (vj,Vi) < Z Vi (vi,Vj) + Z Vj (vi,Vi) + Z Vl (vj,Vi) = ir. LemmaCQthus follows. 

3 The Triangle Witness Algorithm 

In this section, we briefly review the triangle witness algorithm in [6] that constructs the visibility 
graph G = (V, E) of the unknown simple polygon P. 

The triangle witness algorithm is based on Lemma [TJ The algorithm has \n/2\ iterations. In 
the fc-th iteration (1 < k < |7i/2~|), the algorithm checks, for each i = 0, l,...,n — 1, whether Vi 
is visible to Vi + k- After all iterations, the edge set E can be obtained. To this end, the algorithm 
maintains two maps F and B: F[vi\[vj\ = t if Vj is identified as the t-th visible vertex to Vi in 
the CCW order, i.e., Vj = vist{vi)\ the definition of B is the same as F. During the algorithm, 
F will be filled in the CCW order and B will be filled in the clockwise (CW) order. When the 
algorithm finishes, for each Vi, F[vi] will have all visible vertices to Vi on the chain c/i(wj+i, v\ n /2\ ) 
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while B[vi] will have all visible vertices to Vi on the chain ch(yr n /2i ,"Uj_i). Thus, F[vi\ and B[vi] 
together contain all visible vertices of P to V{. For ease of description, we also treat F[v{] and 
B[vj\ as sets, e.g., vi G means that there is an entry and l-Fffi]! means the number 

of entries in the current -F[fj]. 

Initially, when k = 1, since every vertex is visible to its two neighbors along the boundary 
of P, we have = 1 and B\vj\ [wj-i] = deg(vi) for each Uj. In the k-th iteration, we 

determine for each v^, whether Vi is visible to v^k- Below, we let j = i + k. Note that |-F[vj]| + 1 
is the index of the first visible vertex to Vi in the CCW order that is not yet identified; similarly, 
{deg{vj) — \B[vj}\) is the index of the first visible vertex to Vj in the CW order that is not yet 
identified. If Vi is visible to Vj, then we know that Vj is the (|.F[vj]| + l)-th visible vertex to Vi 
and Vi is the (deg(vj) — \B[vj]\)-th visible vertex to vj, and thus we set = + 1 and 

-B[vj][fi] = deg(vj) — \B[vj]\. If Vi is not visible to Vj, then we do nothing. 

It remains to discuss how to determine whether V{ is visible to v j . According to Lemma HJ we 
need to determine whether there exists a triangle witness vertex vi in the chain ch(vi+\, v j-i), i.e., 
vi is visible to both Vi and vj and /^(^z, Vj)+/Xj (vi, vi)+Z Vl (vj, v i) = ir. To this end, the algorithm 
checks every vertex v\ in ch(vi + i, For each vi, the algorithm first determines whether vi is 

visible to both Vi and Vj, by checking whether there is an entry for v\ in F[i>i\ and in The 
algorithm utilizes balanced binary search trees to represent F and B, and thus checking whether 
there is an entry for v\ in F[vj\ and in B[j] can be done in O(logn) time. If vi is visible to both Vi 
and vj, then the next step is to determine whether Zl^v^Vj) + Zj,Jvi,vi) + Z Vl (vj,Vi) = n. It is 
easy to know that Z Vl (vj,Vi) is Z Vl (F[vi] [vj], B[v{\ [vi]), which can be found readily from the input. 
To obtain Zl^v^Vj), we claim that it is the angle /^(F^Jfuj], + Indeed, observe that all 

visible vertices in the chain c/i(t>j+i , Uj-i) are in the current F[vi] . As explained above, \F[vi] | + 1 is 
the index of the first visible vertex to Vi in the CCW order that has not yet been identified, which is 
the first visible vertex to V{ in the chain ch(vj,Vi), i.e, the vertex v'y Thus, Z^(F[^][?;;], |-F[u«]| + 1) 

is Z Vi (vi,v'j), which is zl^v^vj). Similarly, ^(v^vi) is the angle Z Vj (deg(vj) - \B[vj]\, B[vj][vi]). 
Note that both the angles /^(Ffwj]^], |-F[i>i]| + 1) and Z Vj (deg(vj) — \B[vj]\, B[vj][v{\) are known. 
Algorithm [2] in the Appendix summarizes the whole algorithm [6]. 

To analyze the running time of the above triangle witness algorithm, note that it has \n/2] 
iterations. In each iteration, the algorithm checks whether Vi is visible to Vi + k for each < i < n—1. 
For each Vi, the algorithm checks every v\ for % + 1 < / < i + k — 1, i.e, in the chain c/i(ui+i, Uj+fc-i). 
For each such vi, the algorithm takes O(logra) time as it uses balanced binary search trees to 
represent the two maps F and B. In summary, the overall running time of the triangle witness 
algorithm in [6] is 0(n 3 log n). 

4 An Improved Triangle Witness Algorithm 

In this section, we present an improved solution over the triangle witness algorithm in [6] sketched 
in Section Our improved algorithm runs in 0(n 2 ) time. Since the input size (e.g., the total 
number of visibility angles) is 0(ra 2 ) in the worst case, our improved algorithm is worst-case 
optimal. Our new algorithm follows the high-level scheme of the triangle witness algorithm in [6], 
and thus we call it the improved triangle witness algorithm. 

As in [6], the new algorithm also has \n/2] iterations, and in each iteration, we determine 
whether Vi is visible to Vi+k for each < i < n — 1. For every pair of vertices V\ and Uj+fc, let 
j = i + k. To determine whether Uj is visible to Vj, the triangle witness algorithm [6] checks each 
vertex vi in the chain ch(vi+i,Vj—i) to see whether there exists a triangle witness vertex. In our 
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Figure 4: Illustrating Lemma[2] v' is visip! v Ai(vi), the last visible vertex to Vi in ch(vi+i,Vj-i); Vi is visible to Vj 
if and only if v is a triangle witness vertex. 




Figure 5: Illustrating the visibility polygon P(i) for the example in Fig. [3] The gray (pink) area is P(i), which is 
a star-shaped polygon with the vertex Vi as a kernel point. 

new algorithm, instead, we claim that we need to check only one particular vertex, vis\Ft v .-\\(vi), 
i.e., the last visible vertex to V{ in the chain ch(vi + i, v in the CCW order, as stated in the 
following lemma. 

Lemma 2 The vertex v-i is visible to Vj if and only if the vertex vis\p\ v M(vi) is a triangle witness 
vertex for vi and Vj . 

Proof: Let v' denote the vertex vis^pr Vi u(vi), which is the last visible vertex to Vi in the chain 
ch(vi+\, Vj-i) in the CCW order. Recall that v' being a triangle witness vertex for Vi and Vj is 
equivalent to saying that v' is visible to both Vi and vj and Z| 4 (v ', Vj) + Zj. (vi, v') + Z^' (vj, Vi) = it. 

If Vi is visible to Vj, then we prove below that v' is a triangle witness vertex for vi and Vj, i.e., we 
prove that v' is visible to both Vi and Vj and zX^v' ,Vj) + /X j {vi,v') + A v i{vj,Vi) =tt. Refer to Fig. 0] 
for an example. Let P(vi) denote the subpolygon of P that is visible to the vertex V{. Usually, P(vi) 
is called the visibility polygon of i>i and it is well-known that P{vi) is a star-shaped polygon with u, 
as a kernel point (e.g., see pLj). Figure [5] illustrates P(vi). Since v' is the last visible vertex to Vi in 
ch(vi+i,Vj—i) and Vj is visible to Wj, we claim that Vj is also visible to v' . Indeed, since both v' and 
Vj are visible to Vi, v' and Vj are both on the boundary of the visibility polygon P{vi). Let P'ivi) 
denote the portion of the boundary of P{vi) from v' to Vj counterclockwise (see Fig. [6|). We prove 
below that P'[vi) does not contain any vertex of P except v' and Vj. First, P'{vj) cannot contain 
any vertex in the chain ch(y j+i, (otherwise, Vj would not be visible to vi). Let the index of the 
vertex v' be I, i.e., v' = vi. Similarly, P'(vi) cannot contain any vertex in the chain ch(vi+i, 
(otherwise, v' would not be visible to Finally, P'ivi) cannot contain any vertex in the chain 
ch(vi + \,Vj-i), since otherwise v' would not be the last visible vertex to Vi in ch(vi + i,Vj^i). Thus, 
P'ivi) does not contain any vertex of P except v 1 and Vj. Therefore, the region bounded by P'[vi), 
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Figure 6: Illustrating P'(i) for the example in Fig. [4] The thick line segments form P'(i). 



the line segment connecting V{ and Vj, and the line segment connecting Vi and v' must be convex 
(in fact, it is always a triangle) and this region is entirely contained in P. This implies that Vj 
is visible to v' . Hence, the three vertices Vi,Vj, and v' are mutually visible to each other, and we 
have Zl t (vi,Vj) + Zl^vi^i) + Z Vt (vj,Vi) = Z v .(yi, Vj) + Z Vj (v{, v t ) + Z Vl (vj, vi) =tt. 

On the other hand, if v' is a triangle witness vertex for V{ and Vj, then by Lemma [TJ the vertex 
Vi is visible to Vj. The lemma thus follows. □ 

By Lemma [2j to determine whether V{ is visible to Vj , instead of checking every vertex in 
ch(vi + i,Vj-i), we need to consider only the vertex vis\ F [ v .^(vi) in the current set F[vj\. Hence, 
Lemma[2]immediately reduces the running time of the triangle witness algorithm by an 0{n) factor. 
The other O(logn) factor improvement is due to a new way of defining and representing the maps 
F and B, as elaborated below. In the following discussion, let v' be the vertex vis\pi Vi \\(vi) . 

In our new algorithm, to determine whether Vi is visible to Vj, we check whether v' is a triangle 
witness vertex for Vi and Vj. To this end, we already know that v' is visible to Vi, but we still 
need to check whether v' is visible to Vj. In the previous triangle witness algorithm [6], this step 
is performed in O(logn) time by representing F and B using balanced binary search trees. In our 
new algorithm, we handle this step in 0(1) time, by redefining F and B and using a new way to 
represent them. 

We redefine F as follows: i* 1 ^]^] = t if Vj is the i-th visible vertex to Vi in the CCW order; if 
Vj is not visible to Vi or Vj has not yet been identified, then -F[uj][i>j] = 0. For convenience, we let 
F[vi] = F[vi] [uj+i, • • • j Vi-i]- Thus, in our new definition, the size of F[vi\ is fixed throughout 
the algorithm, i.e., I-Fffj]] is always n — 1. In addition, for each Vi, the new algorithm maintains 
two variables Lj and Ij for -F[w«], where Lj is the number of non-zero entries in the current -F[i>i], 
which is also the number of visible vertices to V{ that have been identified (i.e., the number of 
visible vertices to Vi in the chain ch(vi+i, Vi+k—i)) up to the k-th iteration, and I{ is the index 
of the last non-zero entry in the current i.e., Ij is the last visible vertex to Vi in the chain 

ch(vi+i, v i+k—i) m the CCW order. Similarly, we redefine B in the same way as F, i.e., for each 
Vi, B[vi] = B[vi][vi + i,Vi + 2, ■ ■ ■ ,Vi-i) and B [uj] [ty] = t if Vj is the i-th visible vertex to V{ in the 
CCW order. Further, for each Vi, we also maintain two variables L\ and l[ for -B[fj], where L\ is 
the number of non-zero entries in the current -B[uj], which is also the number of visible vertices to 
Vi in the chain ch(vi-k+i, Vi-i) (up to the k-th iteration), and l[ is the index of the first non-zero 
entry in the current -B[uj], i.e., l[ is the first visible vertex to V\ in the chain ch(vi—k+i, i>i-i) in 
the CCW order. During the algorithm, the array F[vi] will be filled in the CCW order, i.e., from 
the first entry F^j^+i] to the end while the array B[vi] will be filled in the CW order, i.e, from 
the last entry Bfujjfttj-i] to the beginning. When the algorithm finishes, F[vi] will contain all the 
visible vertices to Vi in the chain ch(vi+i, Uj+r n /2] )> an d thus only the entries of the first half of 
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Algorithm 1: The Improved Triangle Witness Algorithm 



Input: n, (v±,V2, ■ ■ ■ ,v n ), deg(vi), and the angles Z v< (-, •) for each < i < n — 1. 
Output: The edge set E of the visibility graph G for the target polygon P. 

/* All indices below are understood modulo n. */ 

L^0; 

for i «— to n — 1 do 

F[vi][vi+i,Vi+2, ■ ■ ■ ,Vi-i] <- 0; 
_ B[vi][v i+1 ,Vi + 2, ■ ■ .,Vi-i] <- 0; 

for i <- to n - 1 do 

E ^ EU{e( Vi ,v i+l )}; 
F[vi][v i+1 ] <- 1; 
^^ij^i-i] <r- deg(vi); 
_ Li «- «- Vi+i^L'i <- 1, 1[ <- Vi-i; 

for fc <r- 2 to [n/2] do 
for i to n — 1 do 

j <r- i + fe; 

«' <- 

if B[«j][t/] > then 

4( 1) > J )^4(FMM,L ! + 1); 
Z^y) <- Z Vj (deg( Vj ) - L[, B[vj)[v']); 
\_ Z v ,( Vj ,Vi) <- /VWN.-BMM); 

if Zi< («',«,•) + Zjj.^i/) + Z u /(v 3 -,«i) = 7r then 

F^]^-] <-Li + l; 
B[vj\[vi) <- deg(vj) - L'-; 
|_ Li <- Li + 1, L^ <- L^ + 1, Ji <- /j <- Vi; 



F[vi] are possibly filled with non-zero values. Similarly, only the entries of the second half of B[vi] 
are possibly filled with non-zero values. Below, we discuss the implementation details of our new 
algorithm, which is summarized in Algorithm [TJ 

Initially, when k = 1, for each Vi, we set Ffuijfuj+i] = 1 and £? [^j] = deg(vi), and set all 
other entries of F[vi] and B[vj\ to zero. In addition, we set Lj = 1, L = Vi + \ and L\ = 1, /• = 
In the fc-th iteration, with 1 < k < \n/2], for each m, we check whether Vi is visible to Vj, with 
j = i + k. If Vi is not visible to Vj, then we do nothing. Otherwise, we set L'fuijfvj] = Lj + 1 and 
increase Lj by one; similarly, we set B[uj][«i] = deg(vj) — V- and increase L'- by one. Further, we 
set L = Vj and V- = Uj. 

It remains to show how to check whether m is visible to Vj. By Lemma [21 we need to determine 
whether v' is a triangle witness vertex for and Vj. Since v' is the last visible vertex to Vi in the 
chain ch(vi+i,Vj-i) in the CCW order, based on our definition, v' is the vertex L. After knowing 
v', we then check whether v' is visible to Vj, which can be done by checking whether L^Uj] [i/] 
is zero, in constant time. If jB [?/] is zero, then v' is not visible to vj and v' is not a triangle 
witness vertex; otherwise, v 1 is visible to Vj. (Note that we can also check whether L[f'][uj] is 
zero.) In the following, we assume that v' is visible to Vj. The next step is to determine whether 
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Zl^v'jVj) + ZlJvi,v') + Z v >(vj,Vi) = 7r. To this end, we must know the involved three angles. 
Similar to the discussion in Section E2 we have Z v i(vj,Vi) = Z v >(F[v'][vj], B[v'][vi\), Z.l t (v',Vj) = 
Z v .(F[vi][v'], Li + 1), and = Z Vj (deg(vj) — L'-,B[vj}[v'\). Thus, all these three angles 

can be obtained in constant time. Hence, the step of checking whether V{ is visible to Vj can be 
performed in constant time, which reduces another O(logn) factor from the running time of the 
previous triangle witness algorithm in [6|. Algorithm Q] summarizes the whole algorithm. Clearly, 
the running time of our new algorithm is bounded by 0(n 2 ). 

Theorem 1 Given the visibility angles and an ordered vertex sequence of a simple polygon P, the 
improved triangle witness algorithm can reconstruct P (up to similarity) in 0(n 2 ) time. 

As discussed in [6], the above algorithm can also be used to determine whether the input angle 
data are consistent. Namely, if there is no polygon that can fit the input angle data, then the 
algorithm in Theorem Q] can be used to report this case as well. 
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Appendix 



Algorithm 2: The Triangle Witness Algorithm [6j 

Input: n, (i>i,i>2, • • • , v n ), deg(vi), and the angles Z Vi (-, •) for each < i < n — 1. 
Output: The edge set E of the visibility graph G for the target polygon P. 

/* All indices below are understood modulo n. */ 

F <— [array of n empty maps]; 
B [array of n empty maps]; 

for i <— to n — 1 do 

Ei-EU{e(v i ,v i+1 )}; 
F[vi][v l+1 ] <- 1; 
_ B[vi][vi-i] <- deg(vi); 

for k <- 2 to fn/2] do 
for i <(— to n — 1 do 

j <- i + fe; 

for i •<— i + 1 to j — 1 do 

if vi G and i>; G -B[t>j] then 

4, <- Z Vi (F[vi][vi], \F[ Vi ]\ + 1); 

Z^(«i,uj) <- Z Vj (deg(vj) - \B[vj]\, B[vj)[vi]); 
|_ 4,, <- Z^^F^]^],^^]^]); 

if /l t (vi,Vj) + Z^(wi,u/) + Z Vl (vj,Vi) = it then 
E EU {e(vi,Vj)}; 
F[vi][v3]^\F[vi]\ + l; 
B [vj][vi] <- rfeg(vj) - |Bfe]|; 
abort the innermost loop; 
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